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Abstract. We present a survey of all 3CR sources imaged with ISOCAM onboard the Infrared Space Observatory 
(ISO). The sample consists mostly of radio-loud active galactic nuclei (AGN). For each source, we present spatially 
integrated mid-infrared (MIR, 5 — 18/im) fluxes measured from newly calibrated ISOCAM images. In total, we 
detected 68 objects of the 3CR catalogue, at redshifts z < 2.5, and obtained upper limits for 17 objects. In 
addition, we detected 10 galaxies not listed in the 3CR catalogue. The one with the highest redshift is 4C+72.26 
at z = 3.53. ISOCAM data are combined with other photometric measurements to construct the spectral energy 
distribution (SED) from optical to radio wavelengths. The MIR emission may include synchrotron radiation of the 
AGN, stars of the host galaxy or dust. Extrapolation of radio core fluxes to the MIR show that the synchrotron 
contribution is in most cases negligible. In order to describe dust emission we apply new radiative transfer models. 
In the models the dust is heated by a central source which emits photons up to energies of IkeV. By varying 
three parameters, luminosity, effective size and extinction, we obtain a fit to the SED for our objects. Our models 
contain also dust at large (several kpc) distance from the AGN. Such a cold dust component was neglected in 
previous computations which therefore underestimated the AGN contribution to the far infrared (FIR). In 53 
cases (~ 75 % of our detected 3CR sources), the MIR emission can be attributed to dust. The hot dust component 
is mainly due to small grains and PAHs. The modelling demonstrates that an AGN heating suffices to explain the 
ISO broad band data, starburst activity is not necessary. In the models, a type 1 AGN is represented by a compact 
dust distribution, the dust is therefore very warm and emission of PAHs is weak because of photo-destruction. In 
AGNs of type 2, the dust is relatively colder but PAH bands are strong. 
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1. Introduction 

More than 20 years ago Kotanyi & Ekers (1979) pointed 
out that for radio-loud elliptical galaxies, the orientation 
of the (optical) dust lane appears preferentially perpen- 
dicular to the radio axis. In unification theories dust plays 
an important role. In all of them, a supermassive black- 
hole accretes gas from a disk and, in radio-loud objects, 
a jet is ejected parallel to the rotation axis. The blackhole 
and the disk are surrounded by a torus of interstellar gas 
and dust which may block light from the central region on 
the way to the observer. The diversity in the appearance 
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of activity types is then explained foremost as a result of 
different viewing angles (Barthel 1989). 

The obscuring dust torus is at the heart of unifica- 
tion theories. It reaches towards the blackhole as close as 
dust can survive (T<1500K). As the torus is not directly 
seen on optical images, its dimensions are at most a few 
hundred parsec. Because of the immense bolomctric lumi- 
nosity emitted from the AGN, the dust in the torus, even 
at a distance of 100 pc, must be very warm (> 100 K), and 
further in even hot. Such hot dust will emit preferentially 
in the mid infrared (MIR), at wavelengths covered by the 
ISOCAM filters. 

Dust in radio galaxies is not restricted to the nuclear 
region, but is ubiquitous in the host galaxy. It absorbs 
and scatters the light, but also radiates thermally as it is 
heated cither by the nucleus or by stars. Of course, the 
diffuse and spatially extended dust is much colder than 
the dust near the AGN. Rapid star formation is also a po- 
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tcntial heating source, besides the AGN. A detailed search 
for dust in the nuclear regions of 120 radio sources of the 
3CR catalogue was carried out by de Koff et al. (2000) 
using optical HST images at 0.7[im center wavelength 
and of 0.1" resolution. In one out of three galaxies, they 
found evidence for dust obscuration, with a large variety 
of morphologies such as disks, lanes or filaments. Dust dis- 
tributed smoothly on a large scale of several kpc is difficult 
to detect in absorption maps but reveals itself by infrared 
emission. 

The dust shroud of an AGN is often not transparent 
and to model the emission spectrum one has to compute 
the radiative transfer. Models for optically thick dusty 
nuclei and galaxies have been carried out in various ap- 
proximations (Pier & Krolik 1993, Laor & Draine 1993 
Krugel & Siebenmorgen 1994, Granato & Danese 1994, 
Efstathiou & Rowan-Robinson 1995, Siebenmorgen et al. 
1997, Nenkova et al. 2002, Popescu et al. 2004). Previous 
AGN dust models underestimated the FIR as they did not 
include dust located at large distances from the center and 
thus missed the cold component. To overcome the deficit 
emission in the FIR, these authors add an additional com- 
ponent to the models and attribute it to starburst activ- 
ity (e.g. Rowan-Robinson 2000). In Section 4, we present 
AGN dust models where large scale dust emission is in- 
corporated and where the geometry is radically simplified 
so that they can be described by only three parameters. 

A cold dust component was discovered in 18 sources of 
the 3CR catalogue already by IRAS (Heckman et al. 1994) 
and more recently, in more objects, by ISOPHOT (Fanti 
et al. 2000, Meisenheimer et al. 2001, Van Bemmel et 
al. 2001, Andreani et al. 2002, Spinoglio et al. 2002, Haas 
et al. 2004). These authors detected also a few sources in 
the MIR. The hot dust responsible for it must be close to 
the nucleus and its detection serves as a further test for 
the unified model hypothesis. This motivated us to sur- 
vey the 3CR galaxies in the MIR by means of ISOCAM 
images available from the ISO archive. The sensitivity of 
ISOCAM is 2 — 3 orders of magnitude greater than IRAS 
at 12/^m. 

In Section 2 and 3, we present the images, pertaining 
photometric measurements as well as results. The spectral 
energy distributions (SED) are displayed in Section 4. To 
interprete them, we apply radiative transfer models. A 
discussion of generic SED properties in the infrared for 
different AGN types is given in Section 5. The conclu- 
sions are given in Section 6. In Appendix [S] we remark 
on individual galaxies (previous indications of dust, X-ray 
properties, jets, presence of high excitation lines, source 
morphology, companions). In Appendix 151 we make for 
each galaxy a note on the radiative transfer models of 
Section 4. Details of the ISOCAM observational set up is 
summarised in Tab. IB.il The fully reduced and astromet- 
rically corrected ISOCAM images are shown as overlays 
on optical images in Fig. [Bj 



2. Survey data 

We present the data of our survey, describe which 
ISOCAM fields were extracted from the ISO archive, 
briefly overview the ISOCAM observing modes, and out- 
line data reduction, source identification and photometric 
procedures. 

2.1. Source selection 

The parent catalogue for our sample is the third update of 
the revised third Cambridge catalogue (3CR) by Spinrad 
et al. (1985). It contains 298 sources of which 195 are clas- 
sified as radio galaxies and 53 as radio quasars. Most of 
the objects in the 3CR catalogue are well observed at op- 
tical and near infrared wavelengths. We searched in the 
ISO archive for ISOCAM (Cesarsky et al. 1996) obser- 
vations made with the long wavelength array (LW). All 
ISOCAM images containing, potentially, at least, one 3CR 
source have been extracted. Altogether, we found 146 such 
ISOCAM observations. Often the same source was ob- 
served in multiple filter band passes. On the fully reduced 
ISOCAM images, we determined and improved the astro- 
metric pointing by identifying objects with accurate co- 
ordinates listed in the SIMBAD Astronomical Database 
or the NASA/IPAC Extragalactic Database (NED). The 
images are presented in Fig. El We also made 10 serendip- 
itous detections of galaxies which are not members of the 
3CR catalogue. They are eye-ball detections at the posi- 
tion of NED coordinates. 

2.2. Observations 

Most of the observations were done in the so called mini- 
raster mode where different ISOCAM detector pixels saw 
the same part of the sky (see ISO Handbook for a descrip- 
tion of observing modes, Leech & Pollock, 2001). In ad- 
dition, there are nine staring and one beam switch obser- 
vation. Three observations were performed in the circular 
variable filter (CVF) scan mode and one in the polarimct- 
ric mode. 

Depending on the lens used in the experiment, the 
pixel scale of the 32 x 32 element detector was 1.5", 
3" and 6" and the total field of view 48 x 48, 96 x 96 
and 180 x 180arcsec 2 , respectively. The most frequently 
applied broad band filter was LW10 with a bandwidth 
from 8-15/^m and central wavelength at 12/mi. In addi- 
tion there are observations with filter: LW1 (4-5, 4.5/zm), 
LW2 (5-8.5, 6.7/zm), LW3 (12-18, 14.3/xm), LW7 (8.5-11, 
9.6/im) and CVF scans (5-17/zm) at spectral resolution 
of A/ A A « 50. In a single observing template often more 
than one filter was selected. 

Parameters of the observing templates are detailed in 
Tab. |Rl] target dedicated time (TDT) by which the ob- 
servations are identified (column 1), celestial coordinate 
of the map center (column 2 and 3), observing date (col- 
umn 4), band pass filter (column 5), lens (column 6), num- 
ber of raster points M (column 7) and raster lines N (col- 
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Table 1. ISOCAM photometry of 3CR sources Table 1. - continued. 
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TDT 


Band 


Flux 


RMS 
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0.2 
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14 
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19 
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23 
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4.5 


30 


2 




58703793 


12.0 


25 


1 




27000606 


6.7 


16 


1 




58703801 


14.3 


30 


2 




27000606 


14.3 


13 


1 


3C033.1 


59702607 


12.0 


18 


1 


3C303.1 


52900115 


12.0 


1.9 


0.2 


3C048 


43901804 


14.3 


59 


4 


3C305 


46300408 


4.5 


12 


1 


3f!06i i 


5620141 1 


12.0 


4.3 


0.3 




12300205 


6.7 


7.1 


2.7 


3C66B 


45304902 


4.5 


14 


1 




51400760 


12.0 


21 


1 




43502724 


6.7 


13 


1 




12300205 


14.3 


27 


5 




45304902 


14.3 


3.8 


1.6 


3C305.1 


71702771 


12.0 


1.5 


0.2 


3C071 


63301602 


4.5 


6860 


340 


3C309.1 


60001472 


12.0 


8.2 


0.5 


("Nocinfis") 

111 VJ V7 1-\J\jKJj j 


63301902 

\) *J 'Jul JUiJ 


6.8 


14030 


700 




60092200 


6.7 


4.4 


0.5 


63302202 


14.9 


50830 


2540 


3C319 


54100619 


12.0 


1.3 


0.2 


3C076 1 


46601603 

I V f \ I V / 1 V ' V7 '7 


4.5 


5.1 


1.0 


3C321 


65800208 


6.7 


12 


1 




46601603 


6.7 


3.0 


0.8 




65800208 


12.0 


27 


1 




466(11 60^ 


1 A 


i n 


1 ^ 




65800208 


14.4 


51 


3 


3C079 


61503807 


12.0 


21 


2 


3C324 


30300612 


6.7 


0.1 


0.2 


3f!083 1 


65901 304 

VJ t7 t7 W 1 <J V/ 7t 


4.5 


31 


2 




30300413 


12.0 


1.7 


0.2 




65901032 


6.7 


47 


3 


3C330 


33600323 


6.7 


1.3 


0.2 




65901304 


14.3 


9.7 


1.5 




60001373 


12.0 


1.7 


0.2 


3C084 


61503617 


6.7 


230 


20 




33600323 


14.3 


2.6 


0.5 
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Table 1. - continued. 



Name 


TDT 


Band 


Flux 


RMS 
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22 
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0.8 
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12.0 
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0.2 



umn 8), step between raster points dM (column 9) and 
raster lines dN (column 10), both in arcsec, and the me- 
dian number of exposures taken on each raster position 
(column 11) together with the integration time per read- 
out (column 12). 

2.3. Data reduction 

The data were reduced with the ISOCAM Interactive 
Analysis (CIA, Ott et al. 1996). We used the default data 
reduction steps of CIA: dark current subtraction, initial 
removal of cosmic ray hits (glitches), detector transient 
fitting, exposure coaddition, flat fielding and flux conver- 
sion to astronomical units. Except for staring modes, the 
coadded images at each raster position were sky projected 
and corrected for field distortion. 

The dark current depends on the orbital position of 
the ISO space-craft and the temperature of the ISOCAM 
detector. The applied correction is based on a model de- 
scribed by Roman & Ott (1998). The deglitching is done 
by following the temporal signal variation of a pixel using 
a multi-resolution wavelet transform algorithm (Starck et 
al. 1997). The response of the detector pixels depends on 
the previous observations and there may be long term hys- 
teresis effects for each detector element after changes of 
the photon flux level. We applied the detector flux tran- 
sient fitting method for ISOCAM data (Coulais & Abergel 
2000). After application of the default deglitchcr some 



Table 2. ISOCAM 3cr upper limits of 3CR sources 
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residuals of cosmic ray impacts were still visible in the 
data. Therefore we applied after the detector flux transient 
correction a second cosmic ray rejection method which 
is basically a multi-sigma clipping of the temporal signal 
(Ott et al. 2000). For raster observations, flat field estima- 
tion can be improved by an iterative method exploiting 
the fact that the same sky area is seen by different de- 
tector pixels (Starck et al. 1999). The method works best 
for highly redundant data where a large number of pixels 
see the same sky. In case of low redundancy rasters and 
for staring observations, we flat -fielded according to the 
standard calibration library. 

2.4. Source photometry 

To determine the source flux, we perform a multi-aperture 
photometry. For each aperture centered on the brightest 
pixel of the source we determine the background as the 
mean flux derived in a 4 pixel wide annulus which is put 
2 pixels away from the greatest (~ 18") aperture. In this 
way, we obtain aperture fluxes that flatten with increas- 
ing aperture radius and approach an asymptotic value. 
For point sources the same procedure is repeated on a 
theoretical and normalised point source image (Okumura 
1998) to find the correction factor of the multi-aperture 
analysis. The correction factor is typically 5%. The deter- 
mination if a source is a point source or not is done by 
eye. The procedure will introduce an error of less than the 
correction factor for slightly extended sources misjudged 
to be a point source. In order to derive the statistical flux 
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Table 3. ISOCAM photometry of galaxies not in the 3CR 
catalogue 



Name 


TDT 


Band 


Flux 


RMS 
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uncertainty, for the source aperture we quadratically add 
the RMS image and the la uncertainty of the background 
estimate. The systematic error of ISOCAM photometry 
is typically <5% and color corrections of order 10% (see 
ISOCAM Handbook, Blommaert et al. 2001). 

3. Results 

We find from the detection statistics of the ISOCAM sur- 
vey of 3CR sources that in our sample there is no bias 
against AGN type. We present SEDs from optical to radio 
"wavelengths for individual sources and discuss potential 
contributions to the ISOCAM fluxes (synchrotron, photo- 
spheric and dust emission). 

3.1. Bias and detection statistics 

The primary goal of our survey is to compile a sample of 
radio-loud AGNs which is not subjects to biases related 
to the type and therefore viewing angle of the AGN torus. 
While this seems to be the case for the parent 3CR cata- 
log (e.g. Urry and Pdovani 1995), the selection of objects 
from the ISOCAM survey was beyond our control. We 
shall therefore compare the flux and redshift distribution 
of the galaxies detected by ISOCAM with those of the 
3CR catalogue. 

In total, we detected 68 sources from the 3CR cat- 
alogue and determined 17 upper limits up to redshifts 
z < 2.5. Furthermore, there are 10 serendipitous detec- 
tions of galaxies in the ISOCAM images which are not 



Fig. 1. Top: Histogram of the distribution of the 178 MHz 
radio fluxes in the 3CR catalogue. The filled part repre- 
sents the objects detected by ISOCAM. Bottom: Radio 
fluxes of 3CR sources at 178 MHz plotted as a function of 
redshift. The filled squares are objects of this paper. 



listed in the 3CR catalogue. One example is the radio 
source 4C+72.26 at redshift z = 3.53. The ISOCAM fluxes 
of detected 3CR sources are listed in Tab.Q]and upper lim- 
its of non-detected 3CR sources in Tab. [21 We also give the 
ISOCAM flux for galaxies which are not members of the 
3CR catalogue in Tab. [3] Notations in Tab.^ an( i Tab. 03 
are: catalogue name (column 1), the TDT number used 
as an identifier of the observation is listed in column 2, 
followed by the central wavelength (/im) of the filter band 
pass (column 3), aperture flux (column 4) and la RMS 
flux error (mJy) of the aperture photometry (column 5). 
Tab. [21 has the same structure except that we specify as 
flux the 3a upper limit. Notes to 3CR galaxies detected 
by ISOCAM are given in Appendix lAl 

In Fig. ^ we compare the radio flux and redshift distri- 
bution of 3CR sources detected with ISOCAM to galaxies 
of the parent catalogue. From the top panel of Fig. ^it be- 
comes evident that the fraction of 3CR galaxies detected 
in our survey is a strong function of radio flux, which 
spans two orders of magnitude. For example of the sources 
with an 178 MHz flux <10 Jy we detect only 15%, whereas 
for sources with an 178 MHz flux stronger than 30 Jy 
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Fig. 2. Detection rate of 3CR sources in our survey as a 
function of redshift. The filled part part of the columns 
refers to 3CR sources detected with ISOCAM where the 
MIR emission can be attributed to dust (see SectQJ. The 
shaded part represents all ISOCAM detections of the 3CR 
catalogue and the upper (blank) histogram shows the red- 
shift distribution of 3CR sources sampled by ISOCAM. 



~ 50% of all 3CR sources are detected. Furthermore, the 
galaxies detected with ISOCAM follow, by and large, the 
redshift-flux distribution of the parent catalogue (lower 
box of Fig. nj, at least, for redshifts larger than about 
0.2. However, at z < 0.2, our subsample seems to show a 
preponderance of radio bright objects. 

The fact that low redshift objects are preferentially in- 
cluded raises the question whether this is due to flux lim- 
its of our survey, i.e. whether the ISOCAM detection rate 
strongly drops as a function of redshift. This is investi- 
gated in Fig. where we present the detection rate of our 
survey as a function of redshift. From Fig. [21 one can see 
that the detection rate of local objects (z < 0.1) is close to 
100%. For z > 0.2, and then almost independent of red- 
shift up to z ~ 1, the detection rate is ^75%. Therefore 
our survey is in this redshift range not limited by the sen- 
sitivity of ISOCAM. We conclude that the AGNs of our 
survey present a fair subsample of the 3CR catalogue for 
redshifts larger than 0.2. 

3.2. Spectral energy distributions 

For each detected 3CR source we compiled photomet- 
ric measurements from optical to radio wavelengths. 
Our primary source of information was the NASA/IPAC 
Extragalactic Database (NED). In addition, we availed 
ourselves of the IRAS fluxes (re-processing raw data with 



XSCANPI, sec NASA/IPAC home page), the ISOPHOT 
fluxes (Fanti et al. 2000, Meisenheimer et al. 2001, Van 
Bemmel et al. 2001, Andreani et al. 2002, Spinoglio et 
al. 2002, Haas et al. 2004), mm/submm continuum data 
(Kriigel et al. 1990, Gear et al. 1994, Steppe et al. 1988, 
Steppe et al. 1995, Best et al. 1998, Cimatti et al. 1998, 
Robson et al. 2001, Haas et al. 2004) and 1.4 GHz fluxes 
from the NRAO VLA Sky Survey (Condon et al. 1998). 
The resulting observed SEDs are shown in Fig. [3 together 
with additional information discussed below. 

The MIR emission of our radio galaxies may come from 
synchrotron radiation of the AGNs, stars within the galax- 
ies and dust emission. We assume that bremsstrahlung 
can be neglected. To identify a synchrotron component at 
ISOCAM wavelengths, we extrapolated radio core fluxes 
of the 3CR sources with a power law or, if that was not 
possible, we assumed a spectral slope as derived from 
global radio properties. The latter include the core emis- 
sion (a few arcsecs in size) and the radio lobes which may 
extend several arcmin. The ISOCAM observations have 
a resolution of ~ 10" and are centered on the AGNs. 
Consequently, one has to use radio core fluxes (Giovannini 
et al. 1988; Hardcastle & Worrall 1999) to avoid an over- 
estimate of the synchrotron component. When such core 
fluxes were not available, we adopted the core-to-total ra- 
dio power correlation by Giovannini et al. (2001) to derive 
some estimate. As can be seen from Fig. [31 the synchrotron 
contribution to the MIR emission is generally small. One 
exception is the BL Lac 3C371, which will be neglected in 
the following. 

We assess the significance of the stars by fitting a red- 
shifted T — 4000 K blackbody to the optical and near 
infrared data and extrapolate it to the MIR. The effective 
temperature was chosen deliberately low to maximise the 
effect of the stars at MIR wavelengths. The further inter- 
pretation of this curve depends on the size of the optical 
relative to the ISOCAM image. When they are compara- 
ble, we directly subtracted from the ISOCAM flux the 
extrapolated stellar flux. When the optical observation 
comprised a much larger region, the stellar component was 
considered as the upper limit of a stellar component within 
the ISOCAM aperture. The stellar component is simply 
added to the emission of the galaxy without any radia- 
tive transfer calculations. Indeed, for some of our sources 
the ISOCAM flux is dominated by stars from the host 
galaxy (see, for instance, the SED of 3C270 in Fig. OJ). 
To quantify the likely contribution of the dust emission to 
the ISOCAM fluxes, we present below radiative transfer 
models of dust enshrouded AGNs. 

4. Radiative transfer in dusty AGN 

Here we outline an AGN dust model which is based on 
three free parameters. This simple model seems sufficient 
to fit the dust emission of the observed SED. After some 
motivation we present the model of the central engine 
and of the dust, the parameters and a library of reference 
SEDs. Then we discuss model properties and the fit to the 
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Fig. 3. Observed spectral energy distributions of the sources of the ISOCAM survey: ISOCAM photometry (filled 
squares), NED data (triangles), submillimeter data (pentagons), radio core fluxes (asterisks). Error bars are la and 
upper limits 3er RMS. Solid lines represent the addition of the photospheric component (dashed lines) and the dust 
emission resulting from the AGN models. 
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Fig. 3. - continued - 

data. The computation of the dust emission together with 
the procedure to consider photospheric and synchrotron 
components provides a tool to specify a likelihood of the 
detection of dust in radio - loud AGN. Notes on fits of 
individual sources are given in Appendix iBl 

4.1. Motivation 

In most 3CR objects, the mid and far IR flux cannot arise 
from stars nor from the radio core because an extrapo- 
lation of either component to the infrared fails by orders 
of magnitude. This confirms the widely held view that 
the MIR/FIR emission of 3CR sources is often dominated 
by the dust (Meisenheimer et al. 2001, Van Bemmel et 
al. 2001, Andreani et al. 2002, Freudling et al. 2003, Haas 
et al. 2004). As the sources are not optically thin, one 
has to somehow compute the radiative transfer. In our 
procedure we self-consistently solve the radiative transfer 
in spherical symmetry with a ray-tracing method. The 
heating source of the dust is the central engine, scatter- 
ing by dust is taken into account. Absorption and scat- 
tering efficiencies as well as asymmetry factors of the 
grains are calculated from Mie theory. The size distribu- 
tion which the dust grains display is approximated by 8 
different grain sizes for carbon and silicate material, re- 
spectively. Collisional heating (Contini & Contini 2003) 
is disregarded. The dust model includes a component of 
tiny and transiently heated particles. Further details can 
be found in Siebenmorgen et al. (1992). Observations per- 
taining to the distribution of dust in 3CR galaxies were 
performed by Martel et al. (1999) using the HST. Most 
galaxies could be resolved and dust was found close to 
the central heating source as well as far (few kpc) from 
it. Their images imply a complex morphology of the dust 
obscuration which defies any simple mathematical descrip- 
tion. To be computational tractable, a radical simplifica- 
tion of the geometry is needed. As there are only few ob- 
servational infrared data points for each galaxy, a sophis- 
ticated source structure that would require for its specifi- 
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cation more parameters than available data points is not 
warranted. We therefore consider spherical symmetry. 

This is, of course, not in agreement with the unified 
scheme where the AGN is surrounded by an accreting disk 
or a torus and the phenomenological type is determined 
by the viewing angle of the observer. Such a configura- 
tion is intrinsically, at least, axisymmetric. However, the 
strongest effect brought about by the torus on the emerg- 
ing spectrum occurs, of course, at visual and ultraviolet 
wavelengths. In the infrared, which we try to reproduce, 
the distance of the dust to the source and its mass are 
more important. Therefore, when barring the short wave- 
lengths, a sphere can still conserve key elements of a multi- 
dimensional object. We note also that the angular (not the 
radial) distribution of the dust is in radiative transfer cal- 
culations only of importance when the source is opaque. 

The dust distribution in the model is a spherical aver- 
age over the true distribution. When applied to fitting an 
observed SED, such a model will tend to overestimate the 
near infrared emission of a type 2 AGN (edge-on viewed 
disk) , and underestimate it for a type 1 (pole-on oriented 
disk). Nevertheless, as viewing angles cover a continuous 
range, in a statistical analysis, one should get about the 
right amount of cold and hot dust as well as an estimate 
of their distances from the central heating source. 

Another reason to apply here spherically symmetric 
models stems from the well known difficulty which disks 
have in explaining the 10/im-silicate band. This resonance 
lies in a spectral region covered by ISOCAM. When viewed 
face-on, disks should show the feature in emission, con- 
trary to what is observed. Spherical models, on the other 
hand, are not heavily plagued by this problem as there is 
always some extinction towards the center. As a third ar- 
gument we point out that for one of the best known AGNs, 
in the Circinus galaxy, a spherical model proves as efficient 
in reproducing the observed infrared SED and the wave- 
length dependence of the infrared sizes (Siebenmorgen et 
al. 1997) as do two-dimensional models (Ruiz et al. 2001). 
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The reason lies in the fact that most infrared data are of 
low spatial resolution. 

Previous AGN torus models (e.g. Pier & Krolik 1993, 
Efstathiou & Rowan-Robinson 1995) do not include the 
dust in the host galaxy. Such models therefore miss the 
cold dust altogether. As the apertures of FIR observations 
are large one has to artifically include in these models a 
cold component (Farrah et al. 2003), whereas we incorpo- 
rate it energetically self-consistent. Finally, the dust dis- 
tribution as derived from HST images (de Koff et al. 2000) 
has so far never revealed a perfect torus like structure. 

4.2. The central engine 

The central engine is approximated by a (small) sphere 
with a power law monochromatic luminosity L v oc t/~ - 7 
in the range from 10A to 2/im. The exact lower limit, Ai ow , 
should not be important as the dust becomes transparent 
(Dwek & Smith 1996). The upper limit, A up , is equally un- 
certain and represents a synchrotron break. Fortunately, 
the choice of Ai ow and A up has little influence on the final 
IR spectrum. 

4.3. The dust model 

The dust consists of large carbon and silicate grains with 
radii, a, between 300 and 2400A and a size distribution 
n(a) oc a" 3 5 ; very small graphites (a = 10A); and two 
kinds of PAHs (30 C and 20 H atoms; 252 C and 48 H 
atoms) which are responsible for various infrared emission 
bands (see Siebenmorgen et al. (2001) and Siebenmorgen 
et al. (2004) for details of the dust model). 

The contribution of the different dust populations to 
the AGN spectrum is studied in Fig. 0] For fixed intrinsic 
luminosity of the AGN L = 10 12 L Q , radius R = lkpc 
and dust extinction Ay = 8 mag we show the SED for 
three dust models: i) large grains only; ii) large grains 
and small graphites; iii) all dust components are included 
(large grains, small graphites and PAHs). Of particular 
interest for the interpretation of our ISOCAM observa- 
tions is the MIR. Here dust emission below 6/itm is due 
to small graphites which become very hot (>1000K) after 
absorption of a single UV photon. The emission between 
6 and 9/xm and near 11.3 and 12.8/im is dominated by 
the PAH bands. Large grains with temperatures ~ 300 K 
contribute to the emission at wavelengths ~ 10/im. At 
FIR/submm wavelengths, the emission is entirely due to 
cold large grains (<50K). In summary the entire MIR is 
dominated by the emission of hot dust which may be either 
large grains heated to high temperatures, small graphites 
or PAH. 

4.4. Model library 

In the spectra displayed in Fig. |3J one can usually make 
out three components: optical and near infrared emission, 
probably of stellar origin (<5/mi), the mid and far infrared 




10 

wavelength (p.m) 



Fig. 4. The influence of the dust population on the in- 
frared spectrum of the AGN. Model parameters are: vi- 
sual extinction Ay — 8 mag, radius R = 1 kpc, luminosity 
L = 10 12 - 5 L at distance of 50Mpc. Model with large 
grains only (dashed), with large grains and small graphites 
(dotted), full model with large grains, small graphites and 
PAH (solid line). 

(5/xm <A < 1000/im), and the radio regime (> 1000/im). 
Our fits concern only the second component which we as- 
sume to be due to dust and to originate in the central, 
few kpc region of galactic nucleus. The optical and radio 
emission come from totally separated areas and may thus 
simply be added to the mid and far IR spectrum. 

To keep the models and their subsequent interpreta- 
tion as simple as possible, we consider only three basic 
parameters: the luminosity, the effective linear size of the 
galaxy, and the dust (or gas) mass. Furthermore, we adopt 
a constant dust density. We may then alternatively use as 
parameters: luminosity L, radius R, and visual extinction 

Ay. 

We build a set of reference models 1 by varying L, Ay 
and R in discrete steps over a wide range. In particular, we 
change the luminosity from 10 8 5 to 10 15 25 Lq, the visual 
extinction from 1 to 128 mag and the outer radius of the 
interstellar cloud that surrounds the AGN from 0.125 to 
16 kpc; the inner radius of the dust shell is chosen to be 
close to the evaporation radius of the grains; a typically 
number is 1 pc. 

4.5. Model properties 

In our models, the central engine of the AGN is the sole 
source for heating the dust. Alternatively, the dust could, 
at least partially, be heated by a starburst activity. The 
influence on the spectrum is investigated in Fig. where 
we compare an AGN with a starburst model. In either 

1 The full set of models is available at: 
http : //eso . org/ ~ rsiebenm/ agn_models 
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Fig. 5. Comparison of radiative transfer computations of 
dusty star bursts (solid) and AGNs (dashed). In both mod- 
els we applied a visual extinction of Ay = 32 mag, lumi- 
nosity of L = 10 12 Lq at distance of 50Mpc, and radius 
of R = 1 kpc. 

model, the visual extinction Ay — 32 mag, the luminosity 
L = 10 12 L Q and the cloud radius is R = 1 kpc. In the 
starburst, all stars have a uniform luminosity and temper- 
ature (L* = 20000L Q , T eff = 20000 K) and are surrounded 
by a local dust shell of constant density 10 4 H-atoms/cm 3 
(for details of the radiative transfer in a starburst, see 
Krugel & Siebenmorgen 1994). We see in Fig. OD that in 
the SED of a starburst the PAH features and the submil- 
limeter emission are considerably stronger. As our sim- 
ple AGN model already fit observed spectra of radio-loud 
quasars (Freudling et al. 2003) we simplify our procedure 
and ignore in the following any starburst component. 

Once we have chosen to use AGN models, there are 
three basic parameters in our modelling. The AGN lumi- 
nosity of the models is in most cases well constrained by 
the data and corresponds to the integrated flux. The total 
extinction determines absorption features within the spec- 
trum whereas the outer radius of the dust can be used to 
balance the relative amount of hot and cold dust for a 
given extinction. 

The influence of the extinction, Ay, on the SED is 
shown in Fig. for three different luminosities and for a 
fixed radius of 1 kpc. For the same luminosity and radius 
a high extinction implies strong submillimeter emission, a 
deep silicate absorption feature at 10/im and highly red- 
dened radiation at shorter wavelengths. When the extinc- 
tion is low, on the other hand, the AGN becomes optically 
visible, the silicate feature disappears and the submillime- 
ter flux is weak. The SED peak may be used as a measure 
of the average dust temperature of the bulk material. For 
the same R and L value a lower Ay implies a SED peak 
at shorter wavelengths as compared to models with high 
Ay. On average, the dust is then much warmer. The bulk 
material becomes, of course, also warmer when the AGN 



Fig. 6. AGN model spectra are shown for visual extinc- 
tions of Ay = 1, 2, 4, 8, 16, 32 and 64 mag and luminosities 
of L = 10 10 - 5 ,10 115 and 10 12 5 L Q at a distance of50Mpc. 
The radius is fixed at R = 1 kpc. 
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Fig. 7. The influence of the radius on the AGN model 
spectra. We apply a distance of 50Mpc. Top: For lumi- 
nosity of L = 10 12 L Q and extinction of Ay = 32 mag we 
present models with R = 0.125, 0.25, 0.5, 1, 2, 4, 8 and 
16 kpc. Bottom: For luminosity of L = 10 11 L© and ex- 
tinction of Ay — 8 mag we show models with R = 0.125 
and 16 kpc. 

luminosity increases: models with same extinction and ra- 
dius show a SED peak shifted to shorter wavelengths at 
higher luminosities (Fig. |BJ). 

In Fig. extinction and luminosity are held constant 
and the influence of the radius, i?, on the SED is illus- 
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Table 4. Infrared dust luminosity of 3CR radio galaxies. Galaxies where spectral shape of the dust emission is 
constrained are marked in bold. 



Name 


z 


Type 


log(L dust ) 


R 


A v 


Dust ? 








[Lo] 


[pc] 


[mag] 




3C006.1 


0.840 


NLRG/FRII 


11.75 


16000 


32 


certain 


3C013 


1.351 


NLRG/FRII 


12.25 


16000 


2 


likely 


3C017 


0.220 


BLRG 


11.5 


125 


128 


likely 


3C018 


0.188 


BLRG 


12 


16000 


32 


certain 


3C020 


0.174 


NLRG/FRII 


11.25 


2000 


16 


certain 


3C022 


0.937 


BLRG/FRII 


12.5 


500 


8 


possible 


3C031 


0.017 


LERG/FRI 


10 


125 


2 


certain 


3C033.1 


0.181 


BLRG/FRII 


11.25 


125 


16 


certain 


3C048 


0.367 


QSO 


12.75 


4000 


32 


certain 


3C061.1 


0.186 


NLRG/FRII 


11 


1000 


16 


certain 


3C66B 


0.212 


LERG/FRI 


9.75 


16000 


2 


certain 


3C071 


0.004 


NLRG 


11.5 


2000 


16 


certain 


3C076.1 


0.032 


LERG/FRI 


<9.75 


16000 


2 


no evidence 


3C079 


0.256 


NLRG/FRII 


11.75 


500 


16 


certain 


3C084 


0.018 


NLRG/FRI 


11.25 


2000 


4 


certain 


3C098 


0.031 


NLRG/FRII 


10 


125 


32 


certain 


3C231 


0.007 


LERG/FRI 


10.5 


1000 


16 


certain 


3C249.1 


0.312 


QSO 


11.75 


125 


32 


certain 


3C265 


0.811 


NLRG/FRII 


12.5 


125 


64 


certain 


3C270 


0.007 


NLRG/FRI 


< 8.75 


2000 


32 


no evidence 


3C272.1 


0.003 


LERG/FRI 


9 


500 


2 


certain 


3C273 


0.158 


QSO 


12.75 


500 


1 


possible 


3C274 


0.001 


NLRG/FRI 


< 9 


250 


2 


no evidence 


3C277.3 


0.086 


BLRG 


10.5 


250 


1 


certain 


3C286 


0.849 


QSO/ 


12.5 


8000 


64 


certain 


3C288.1 


0.961 


QSO 


12.5 


1000 


128 


likely 


3C293 


0.045 


LERG/FRI 


10.75 


8000 


2 


certain 


3C295 


0.461 


NLRG/FRII 


12.25 


8000 


1 


certain 


3C296 


0.461 


LERG/FRI 


10 


16000 


8 


likely 


3C303.1 


0.267 


NLRG/FRII 


11.5 


16000 


1 


certain 


3C305 


0.041 


NLRG/FRI 


10.75 


8000 


1 


certain 


3C305.1 


1.132 


NLRG 


< 12 


8000 


2 


no evidence 


3C309.1 


0.905 


QSO 


12.75 


250 


32 


certain 


3C319 


0.192 


LERG/FRII 


11.25 


8000 


128 


likely 


3C321 


0.096 


NLRG/FRII 


11.5 


2000 


16 


certain 


3C324 


1.206 


NLRG/FRII 


12.25 


4000 


8 


certain 


3C330 


0.550 


NLRG/FRII 


11.75 


125 


32 


certain 


3C332 


0.152 


BLRG/FRII 


11.25 


16000 


2 


certain 


3C336 


0.927 


QSO 


12 


250 


4 


certain 


3C338 


0.029 


NLRG/FRI 


10.25 


16000 


16 


certain 


3C341 


0.448 


NLRG/FRII 


11.75 


125 


2 


certain 


3C343 


0.988 


QSO 


11.75 


2000 


16 


possible 


3C345 


0.593 


QSO 


13 


16000 


64 


certain 


3C346 


0.161 


NLRG/FRII 


11 


16000 


16 


likely 


3C349 


0.205 


NLRG/FRII 


11.25 


16000 


32 


certain 


3C351 


0.372 


QSO 


12.25 


500 


32 


certain 


3C356 


1.079 


NLRG/FRII 


12.25 


4000 


64 


certain 


3C371 


0.051 


BL Lac 


< 11.5 


16000 


32 


no evidence 


3C379.1 


0.256 


RG 


11.5 


4000 


64 


certain 


3C380 


0.691 


QSO 


12.5 


1000 


32 


certain 


3C381 


0.159 


BLRG/FRII 


11.25 


250 


32 


certain 


3C382 


0.058 


BLRG/FRII 


11 


125 


2 


certain 
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Table 4. continued 



Name 




z 


Type 


log(Ldust) 

[Le] 


R 
[pel 


A Y 
[mag] 


Dust ? 


3C386 





.017 


LERG/FRI 


9 


125 


16 


possible 


3C388 





.090 


LERG/FRII 


<10.75 


8000 


64 


no evidence 


3C390.3 





.056 


BLRG/FRII 


11 


125 


4 


certain 


3C402 





.023 


NLRG 


10 


16000 


128 


certain 


3C410 





.249 


RG 


12.75 


16000 


128 


certain 


3C418 


1 


.686 


QSO 


< 8.5 


16000 


128 


no evidence 


3C430 





.054 


NLRG 


< 9 


250 


128 


no evidence 


3C442 





.026 


LERG/FRI 


< 8.5 


16000 


128 


no evidence 


3C445 





.056 


BLRG/FRII 


11.5 


125 


16 


certain 


3C449 





.017 


BLRG 


9.75 


8000 


64 


possible 


3C452 





.081 


NLRG/FRII 


11 


250 


64 


certain 


3C456 





.233 


LERG 


11.5 


250 


64 


certain 


3C459 





.200 


NLRG/FRII 


12.25 


2000 


128 


certain 


3C465 





.030 


LERG/FRI 


10.25 


125 


32 


certain 


3C469.1 


1 


.336 


LERG/FRII 


< 8.5 


16000 


128 


no evidence 


4C+72.26 


3 


.532 


RG 


13.5 


250 


64 


no evidence 



trated. We present models with high extinction Ay = 
32 mag and luminosity L = 10 12 L© as well as moder- 
ate extinction Ay = 8 mag and luminosity L = 10 11 L©, 
respectively. When the radius is large, most of the dust is 
farther away from the AGN and therefore on average much 
cooler when compared to models with small R. When the 
radius of the model is increased while the extinction is 
kept fixed more dust is needed for the same submillimeter 
flux. As the luminosity is constant the MIR emission be- 
comes fainter for models with large R and the SED peak 
is shifted to longer wavelengths. 

Another important feature of the SED is the strength 
of the PAH emission. In Fig. (bottom) one notices that 
PAH band emission increases for models with increasing 
cold dust, respectively large radii. If R is small the models 
predict more hot dust and relatively weak or fully absent 
PAH bands. This is because most of the dust emission 
originates nearby the central source and only few PAH 
survive the harsh AGN environment. 

To constrain the models one needs MIR and FIR data. 
The FIR observations shall give a measure of the SED 
peak emission. We have studied the degeneration of our 
AGN models. As may be obvious, for sources where the 
FIR peak is reasonably well known but there is only a 
single data point in the MIR it is still possible to fit the 
SED by a combination of low extinction and large radius 
or vice versa. However, when more MIR data are avail- 
able the degeneration almost disappears. In particular for 
galaxies where the ISOCAM data are suited to indicate 
the depth of the silicate band the models are quite unique. 
This is because the depth of the silicate absorption is di- 
rectly related to the visual extinction. 



4.6. Fit to Data 

We selected models from the AGN library (see Sect. 14.40 
which are consistent with the data. First we require that 
the models smoothed to the ISOCAM bandwidths pre- 
dict the measured ISOCAM fluxes to within 10% accu- 
racy. Second, we rejected models which violate any 3<r 
upper limit. Of the remaining models, we choose the ones 
which best fit the optical to FIR/submillimeter data in 
a least square sense. The optical data can be reasonably 
approximated by a blackbody of a few thousand K. We 
therefore assume that they are due to (weakly reddened) 
stars and this stellar component is fit separately. As a re- 
sult one obtains the stellar luminosity (see Fig. 8). In 16 
out of 69 cases the Ay of the model is <2mag and then 
the AGN may also contribute at optical wavelengths. The 
model flux (AGN) is then subtracted from the measured 
flux and only the rest is used to evaluate the luminosity 
of the stars. The models are shown as solid lines in Fig. El 
They fit generally well, exceptions are: 3C273, or 3C382 
(see Appendix iBl for notes on individual fits). 

The results of our modelling procedure are summarised 
in Tab. 0| We specify the name, redshift (as available 
in NED), type (Veron-Cetty & Veron 2000, Spinrad et 
al. 1985), infrared dust luminosity (Sect.^J and specify a 
likelihood of dust detection. We distinguish four cases: 

— No evidence for dust emission is found. ISOCAM and 
other flux measurements are consistent with a photo- 
spheric or synchrotron component. 

— Detection of dust emission is possible. ISOCAM pho- 
tometry is only slightly above the photospheric com- 
ponent and there are only far infrared data close to 
their detection or confusion limit. 

— Detection of dust emission is likely. Here ISOCAM 
data are well above the predicted photospheric compo- 
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nent but there is a possible synchrotron contribution. 
There are no firm far infrared detections. 
— Detection of dust emission is certain. There are 
reliable infrared detections well above the predicted 
photospheric component and the contribution by 
synchrotron or free-free emission can be ruled out. 



5. Discussion 

5.1. Energy budget: nucleus versus host 

First, we compare how the AGN and the stars of the host 
galaxy contribute to the dust heating. The influence of a 
central starburst is not considered. One of the three basic 
parameters of each SED model is the AGN luminosity. In 
our spherical models, which are all optically thick (Ay > 
1), the AGN luminosity is entirely re-radiated by dust in 
the infrared and it is denoted Ldust- The upper panel of 
Fig. [8] shows the AGN luminosity of all galaxies detected 
by ISOCAM as a function of redshift. In the lower panel, 
we divided Ldust by the luminosity of all the stars, L sta ,r, in 
the galaxy. The latter was estimated from the black-body 
fits of the photospheric component displayed in Fig.|3and 
are not de-reddened. 

One notes a trend of Ldust with z which is typical for 
flux limited samples. As demonstrated in the lower panel 
of Fig. |SJ the AGN luminosity is much higher than the 
stellar luminosity, L s tar, which is concentrated in the near 
infrared and optical region. The powerful AGN heats the 
dust not only in its vicinity but throughout the galaxy 
and acceptable model fits to the observed SEDs were pos- 
sible only when dust at large (several kpc) distances was 
included. 
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5.2. Dust emission and AGN unification 

For 27 AGNs of our sample, marked bold face in Tab. 
there exist enough photometric data to constrain the spec- 
tral shape of the dust emission. For a more detailed com- 
parison, we distinguish between galaxies of optical spectral 
type 1 and 2 and between objects of high (L > 10 11 L Q ) 
and low (L < 10 11 L Q ) luminosity L. Altogether we thus 
have four classes: 6 BLRGs (type 1), 5 QSOs (type 1), 6 
low L and 10 high L NLRGs (both of type 2). According 
to the unification scheme, type 1 AGNs are viewed nearly 
pole-on and have visible broad line region, whereas type 
2 AGNs are viewed edge-on, so that the large extinction 
from the AGN torus obscures the broad-line region. 

In FigHH we plot for the 27 objects the best fit pa- 
rameters: radius and visual extinction. AGNs of type 1 
occupy the left and those of type 2 the right part of the 
frame, so the models for them differ foremost by the size 
parameter or radius R. This was already evident from the 
infrared SEDs in Fig. [7] (bottom) where small model radii 
correspond to BLRGs and large radii to low luminosity 
NLRGs. For low luminosities, both NLRGs and BLRGs 
are fit by moderate extinction (^4y ~ 7 mag). For high lu- 



Fig. 8. Top: The intrinsic AGN luminosity of the model, 
Ldust, as function of redshift. Solid symbols are type 1 
AGNs: triangles QSOs and squares BLRGs. Open squares 
are NLRGs and upper limits as of Tab. Bottom: 
Luminosity ratio of AGN and stellar components as a 
function of redshift. Symbols as in top panel. 



minosity NLRGs and QSOs, we find higher values of Ay, 
the averages are Ay = 20 and 30 mag, respectively. This 
trend is reflected for QSOs and high luminosity NLRGs 
in the top panel of Fig. where SEDs with a radius close 
to the mean value of both AGN types can be compared. 

In summary, the model spectra of BLRGs and QSOs 
(Fig. [7| bottom) peak at ~ 40/im, whereas high and low 
luminosity NLRGs attain their maxima much longer wave- 
lengths, around 150 /im. Therefore, the dust in AGNs of 
type 1 is warmer than in those of type 2. This difference 
is qualitatively consistent with the unification schemes. 
For type 1 AGNs, which are viewed close to pole-on, the 
broad-line region is visible and unobscured hot dust dom- 
inates the SED. By contrast, towards NLRG, which are 
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Fig. 9. The best fit model parameters, radius and dust 
extinction, for 27 galaxies with well determined SEDs. 
The objects are grouped according to their luminosity 
(high or low) and spectral type (1 or 2). There are QSOs 
(type f, high L, filled triangles), BLRGs (type 1, low L, 
filled squares), NLRGs (type 2, high L, open triangles) and 
NLRGs (type 2, low L, open squares). For each type, log- 
arithmic means are marked by crosses. 



viewed face-on, large extinction from the torus obscures 
the broad-line region and pre-dominantly cooler dust is 
detected. The type 1 AGNs also show much less PAH 
emission than type 2. The weakness of the band emission 
in type 1, compared to type 2, is probably due to PAH 
destruction. Averages of observed SED are presented by 
Frcudling et al. (2003) confirming the trends indicated in 
the models presented here. 

6. Conclusions 

ISOCAM detected 80% (71 out of 88) 3CR sources. We 
find evidence of hot dust in 53 radio galaxies or 75% 
of the sources detected in our survey. For each detected 
source, we compiled photometry from optical to radio 
wavelengths. 

The extrapolation of radio core fluxes to ISOCAM 
wavelengths shows that synchrotron emission from the 
core is usually negligible in the MIR, exceptions being 
flat spectrum radio sources. The contribution of the host 
galaxy to the MIR is also generally small, here excep- 
tions are a few sources of type FRI. Thus for most objects 
the origin of the MIR can be attributed mainly to dust. 
This dust is hot and heated by the central engine. The 
emission comes from large grains with temperatures of a 
few hundred Kelvin and small grains and PAHs under- 
going temperature fluctuations. This picture is supported 
by simple radiative transfer calculations. Most SEDs can 
be successfully fit between 1 and lOOO^rn by a three pa- 
rameter model where AGN luminosity, size of the galaxy 
and extinction are varied. Our modelling shows that the 



broad band data are consistent with AGN only heating of 
the dust. 

For similar luminosities, we compare AGNs of type 1 
and 2. We find that the model parameters depend strongly 
on the AGN type and a dichotomy of the infrared SEDs is 
derived. The IR fluxes of BLRGs and QSOs peak around 
40/im while the maximum of the dust emission for NLRGs 
is reached at >100^m. Therefore, the dust in type 1 AGNs 
is warmer than in type 2. The models also predict much 
weaker PAH emission in type 1 AGNs than in type 2, as 
a result of evaporation. Unfortunately, PAH bands cannot 
be resolved by the broad-band ISOCAM observations. 

AGNs and starbursts are accompanied by tremendous 
IR luminosities. Although the possible coexistence of both 
energy sources and their relative contributions to the dust 
heating is difficult to assess, starbursts tend to favour PAH 
emission and this may serve as a possible discriminator of 
activity. 
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Appendix A: Notes on individual galaxies 

Our remarks below concern the following points: (1) previous 
indications of dust, (2) X-ray emission, (3) jets, (4) high exci- 
tation line emission, (5) morphology, (6) companions. 
3C002 (z — 1.037): Fuzzy optical structure similar to radio 
morphology (Ridgway & Stockton 1997). 

3C006.1 (z = 0.804): There is an optical extension with posi- 
tion angle roughly equal to that of the radio axis (McCarthy 
et al. 1997). 

3C013 (z — 1.351): There are two faint tail-like structures to 
the south and east of the nucleus. The core radio source is a 
double with a size of 28". 

3C017 (z = 0.2): X-ray source (Siebert et al. 1996) with a 
strong and broad Ha line (Eracleous & Halpern 1994). 
3C018 (z — 0.19): Shows a rich emission line spectrum with 
extremely high ionization states (Tadhunter et al. 1993). 
3C020 (z — 0.17): Two cones of emission and a possibly dust 
disk are detected by de Koff et al. (2000) . The large-scale FR II 
class radio structure is aligned roughly perpendicular to the 
proposed dust lane. 

3C022 (z = 0.935): From JKL' images. Simpson et al. (1999) 
conclude that 80% of the light in a 3" aperture is from a weakly 
reddened (Av — 2 — 5mag) quasar; extended near IR light is 
presumably photospheric. 

3C031 (z = 0.0167): Associated with NGC383. HST images 
show a face-on dust disk with a diameter of 7" and an unre- 
solved core (Xu et al. 2000). Martel et al. (1999) confirm an 
optical synchrotron jet. The molecular gas content is M(H2) ~ 
10 9 M (Leon et al. 2001). 

3C33.1 (z = 0.181): Broad line radio galaxy in a region of high 
foreground stellar density (McCarthy et al. 1995). A star is 
superimposed southeast of the nucleus. 



3C04-8 (z — 0.367): Significant part of the OIII line emission is 
confined to the 0.5" central region (Axon et al. 2000). 
3C61.1 (z = 0.186): Narrow-line galaxy. 

3C66B (z = 0.0212): HST observations confirm an optical syn- 
chrotron jet (Martel et al. 1999). Orbital motion in this radio 
galaxy provides evidence for a supermassive black hole binary 
(Sudou et al 2003). 

3C071 (NGC1068, z = 0.003): Famous Seyfert 2; MIR images 
of high spatial resolution are obtained by Siebenmorgen et al. 
(2004). 

3C076.1 (z = 0.0325): A galaxy with no optical or near IR 
features (Colbert et al. 2001). 

3C079 (z = 0.256): Narrow Ha and H/3 lines (Laing et 
al. 1994). 

3C083.1 (NGC1265, z = 0.025): Edge-on radio galaxy with 
radio jet and nearly orthogonal narrow dust lane of ~ 1.5 kpc. 
There is a bright foreground star near the center (Martel et 
al. 1999). 

3C084 (NGC1275, z = 0.0176): Extensively studied giant el- 
liptical galaxy. It has filamentary dust structures which extend 
out to 17 kpc (de Koff et al. 2000). 

3C098 (z = 0.031): Seyfert 2 with extended X-ray emission 
and a compact central X-ray source (Hardcastle & Worrall 
1999). No optical/near IR features (Colbert et al. 2001). 
3C231 (M82, z = 0.001): Archetype starburst galaxy; see 
Kriigel & Siebenmorgen (1994) for a model of the dust emis- 
sion. 

3C249.1 (PG 1100+772, z = 0.312): Absorption line quasar. 
3C264 (z — 0.022): An apparent optical "ring" is claimed to 
be produced by a circumnuclear dust disk (Baum et al. 1997). 
3C265 (z = 0.811): Complex optical structure without a point 
source at the nucleus. Probably highly obscured AGN (Av > 
15mag, Simpson et al. 1999). 

3C270 (NGC4261, z = 0.007): Quasar with nuclear dust 
disk perpendicular to the large-scale radio jets (Jaffe et 
al. 1996). Indication of vigorous star formation around the 
AGN (Canalizo & Stockton 2000). The quasar suffers little 
extinction (Leon et al. 2001). 

3C272.1 (M84, z = 0.004): Elliptical galaxy with dust extend- 
ing out to almost 1 kpc (de Koff et al. 2000). 
3C273 (z = 0.158): First discovered QSO (Schmidt 1963). 
3C274 (M87, NGC4486, z = 0.004): It has an ionized disk 
surrounding the central black hole (Ford et al. 1994, Harms et 
al. (1994) of 2.6 • 10 9 M Q (Lauer et al. 1992). 
3C277.3 (Coma A, z = 0.086): Elliptical with radio jet. 
3C286 (z = 0.849): Contains high excitation line emission 
(Boisse et al. 1998). 

3C288.1 (z = 0.961): Double radio quasar (Akujor et al. 1994). 
3C293 (UGC 08782, z = 0.045): Has a double nuclei core with 
a steep radio spectrum and filamentary dust lanes (Martel et 
al. 1999). 

3C295 (z = 0.464): Extended double-lobe radio galaxy without 
a core (Akujor et al. 1994, Fomalont et al. 2000). 
3C296 (NGC5532, z = 0.0237): This elliptical galaxy hosts a 
sub-kiloparsec, elongated and sharp dusty structure centerd on 
its nucleus (Martel et al. 1999). 

3C303.1 {z = 0.267): The nucleus has an elongated shape. The 
extended radio source is an asymmetric double ~ 1.8" in size, 
elongated perpendicular to two dust lanes (de Koff et al. 2000). 
3C305 (z — 0.039): An inclined disc of gas and dust encircling 
a quasar reddened by Ay > 4 has been observed by Jackson 
et al. (2003). 
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3C305.1 (z = 1.132): The NED coordinates have an uncer- 
tainty of 4". In the DSS, a contaminating source is visible 8" 
away. The radio emission of this high-redshift galaxy forms a 
core-free structure with two lobes separated by ~ 2". 
3C309.1 ([HB89] 1458+718, z = 0.905): This quasar has a 
compact steep radio spectrum (Peng et al. 2000) and remains 
unresolved in the optical (de Vries et al. 1997). 
3C319 (z = 0.192): Resembles a normal non-star-forming el- 
liptical with an optical long axis almost perpendicular to the 
radio axis (Roche & Eales 2000). 

3C321 (z = 0.096): Blue, high surface brightness galaxy. The 
disturbed dust lane with a projected length of ~ 3.7" breaks 
in two components belonging to two merging galaxies (Roche 
& Eales 2000). 

3C324 (z = 1.206): A gravitationally-lensed, extremely elon- 
gated high redshift radio galaxy (Chambers et al. 1996) with 
a steep radio spectrum (a ~ 1.07, Brunetti et al. 1997). It is 
a ROSAT source with an X-ray luminosity of 4.4 x 10 10 L Q 
(Hardcastle & Worrall 1999). From SCUBA and Effelsberg ob- 
servations Best et al. (1998) conclude that the submillimeter 
signal is not associated with synchrotron emission. 
3C330 (z = 0.550): Optically very luminous and relatively 
blue galaxy. The disk structure is dominated by line emission 
(McCarthy et al. 1996). 

3C332 (z = 0.152): This broad line radio galaxy shows in HST 
images a smooth elliptical nuclear region (Taylor et al. 1996). 
3C336 (PKS 1622+23, z = 0.927): This source has a large 
angular radio size (28"). The brightest ISOCAM pixel is close 
to the radio source FIRST J162438.8+234505. The latter may 
be a resolved component of PKS 1622+23. 
3C338 (NGC6166, z = 0.030): Elliptical in the cluster 
Abell 2199. An arc-like dust feature extends for ~ 3" from 
the nucleus to the West (Capetti et al. 2000) . 
3C341 {z = 0.448): As the NED coordinates are uncertain by 
3", we use coordinates derived from DSS2. Martel et al. (1999) 
confirm an optical synchrotron jet. 

3C343 (z = 0.988): An optically elongated and not very cen- 
trally concentrated quasar. The optical and radio emissions 
may have similar shapes (de Vries et al. 1997). 
3C343.1 (z = 0.75): No radio core has been detected (Ludke 
et al. 1998). 

3C345 (z = 0.595): A supcrluminal optically violently vari- 
able quasar (Penston & Cannon 1970, Katajainen et al. 2000). 
The spectrum is flat for wavelengths > 3mm and steepens for 
A < 3 mm (Rantakyroe et al. 1998). As NED coordinates are 
uncertain by 3" we use for astrometry the VLA map by Leahy 
& Perley (1991). 

3C346 (z = 0.162): Narrow-line radio galaxy showing both 
extended and compact X-ray emission (Hardcastle & Worrall 
1999). In the optical there is a jet and a bright core (de Koff et 
al. 1996). We ignore deviant fluxes by Morabito et al. (1982). 
3C349 (z — 0.205): There are two companions, a faint filament 
suggests an interaction within the triple (Roche & Eales 2000). 
3C351 (z = 0.372): Radio- loud QSO of moderate X-ray lumi- 
nosity. The spectrum from 0.3-10. l/u,m can be fitted with a 
broken power law (Neugebauer et al. 1987). 
3C356 (z — 1.079): In the optical, there are two nuclei 5" apart 
(McCarthy et al. 1995, 1997) and both are detected in the 
radio (Eales & Rawlings 1990, Fernini et al. 1993). We identify 
the one which is brighter in the K band with the nucleus of 
the radio source. 3C356 is detected at X-rays with ROSAT 
(Hardcastle & Worrall 1999, Crawford & Fabian 1996). 



3C368 (z — 1.131): The spectrum from radio to optical wave- 
lengths can be fit by a power-law with a break in the infrared. 
There are two knots emitting synchrotron emission in the opti- 
cal (O'Dea et al. 1999). Near infrared spectroscopy by Jackson 
& Rawlings (1997) shows that the nucleus is contaminated by 
a galactic M star. 

3C371 (z = 0.051): BL Lac object with an unresolved HST 
point source (Martel et al. 1999). 

3C379.1 (z = 0.256): Shows extended OBI emission over a 
scale of 9" (McCarthy et al. 1995). 

3C380 (z — 0.692): There are two optical hot spots coinciding 
with radio peaks (Ludke et al. 1998). Therefore, the optical hot 
spots probably originate from synchrotron radiation (de Vries 
et al. 1997). 

3C381 (z — 0.161): Shows extended OIII emission on a scale of 
~ 20" (McCarthy et al. 1995). There is a low surface brightness 
companion 5.5" away (Roche & Eales 2000). 
3C382 (z — 0.058): A classic broad- lined radio galaxy whose 
optical luminosity is dominated by an extremely blue AGN 
(Roche & Eales 2000). 

3C386 (z — 0.017): An elliptical with a compact optical nu- 
cleus. There are two point sources of similar magnitude ~ 3" 
away from the nucleus. They might be faint foreground stars 
or small companion galaxies (Martel et al. 1999). 
3C388 (z = 0.091): This radio galaxy has an unresolved optical 
nucleus (Martel et al. 1999) and two companions with a nuclear 
separation of 7" (Roche & Eales 2000). The X-ray images show 
significant diffuse emission (Hardcastle & Worrall 1999). 
3C390.3 (z = 0.056): A double-lobed, broad line blazar 
(Fomalont et al. 2000) detected also at X-rays (Hardcastle & 
Worrall 1999). The optical band is dominated by a nuclear 
point source (Martel et al. 1999). 

3C401 (z — 0.201): A radio galaxy with a double nucleus (sep- 
aration 3.6") and several neighbours (Roche & Eales 2000). 
3C402 (z = 0.025): An elliptical galaxy with an FR I morphol- 
ogy and a double radio nucleus. There is a massive nuclear disk 
of gas and dust (De Juan et al. 1996). 

3C410 (z = 0.249): An elliptical with an optical jet candidate 
(de Koff et al. 1996). 

3C418 (z = 1.69): Has several nearby but faint companions 
(Lehnert et al. 1999). 

3C430 (z = 0.056): An elliptical with a narrow, arc-shaped 
dust lane that slices the core (Martel et al. 1999). 
3C433 (z = 0.101): For the near infrared we use the data of 
2MASXi J2123445+250427. The absorption map by de Koff et 
al. (2000) shows patchy dust features. 

3C442 {z = 0.026): A galaxy pair interacting with NGC7237 
(De Vaucouleurs et al. 1976). 

3C445 (z — 0.056): A broad emission line object (Eracleous & 
Halpern 1994) with an unresolved nucleus (Martel et al. 1999). 
3C449 (z = 0.017): A low-power radio galaxy with multi- 
ple X-ray components and a nearby companion (Worrall & 
Birkinshaw 2000). Leon et al. (2001) derived from CO a molec- 
ular gas mass of M(H 2 ) ~ 2.3x10 s M Q . The HST images re- 
veals a 5" inner disk-like feature (Zirbel & Baum 1998). 
3C452 (z = 0.081): dc Koff et al. (2000) find indication of a 
faint dust lane near the nucleus. The radio source has a bright 
core and two lobes. 

3C454-1 (z = 1.841): This source appears in the radio as a 
double with no core (Ludke et al. 1998). In the optical con- 
tinuum one sees two nuclear condensations separated by 2" 
(Djorgovski et al. 1988). 
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3C456 (z = 0.233): The radio core of this FR II galaxy lies 
within 1" of the optical core position (Harvanek & Hardcastle 
1998). The near IR shows extended emission (de Vries et 
al. 1997). 

3C459 (z = 0.02): This steep spectrum ULIRG (IRAS 
23140+0348) shows optical absorption features suggestive of 
young stellar populations (Miller 1981) and broad wings on 
both Ha and H/3 (Yee & Oke 1978). In the optical there are 
two polarisation components (Draper et al. 1993) none of which 
coincide with the radio axis. The radio emission (Condon et al. 

1998) is extended to only ~ 10" (48kpc). 

3C465 (NGC7720, z = 0.030): HST reveals a dusty kpc-size 
disk (Martel et al. 1999) which is roughly perpendicular to the 
radio jets (de Koff et al. 2000). 

3C469.1 (z = 1.336): There is a double nucleus in the optical 
not aligned with the radio axis (McCarthy et al. 1997) . 
4C+72.26 (z = 3.536): The galaxy with the greatest red shift 
in our sample. 

NGC5532NED02 (z ~ 0.0237): As no redshift is found in NED 
we use the value for NGC5532 (3C296). 

NGC7236 (z = 0.026): A peculiar radio galaxy merging with 
NGC7237. It is also detected at X-ray (Hardcastle & Worrall 

1999) . 

PMN J0214-1158 (z = 2.34): A large high redshift radio source 
with a steep radio spectrum; simple double morphology and jet 
signatures (Roettgering et al. 1994). 

IRAS F 17130+5021 (z = 0.05): A galaxy pair 15" away asso- 
ciated with 2MASXi J 1714167+501816. 

Appendix B: Notes on individual fits 

3C006.1: Evidence for dust comes from ISOCAM 12^m (this 
work) and weak millimeter detection by Haas et al. (2001). 
Evidence for hot dust comes from ISOCAM only. 
3C013 and 3C017: Fitting the visible and ISOCAM data points 
by a Planck function requires a rest frame temperature of 
2600 K. If one allows for moderate reddening of Ay = 1 mag 
the temperature would rise to 2900 K. Such a value is still lower 
than the expected mean stellar temperature in the galaxy and 
therefore the presence of hot dust is likely. 
3 CO 18: 

3C018: IRAS fluxes suggests Av larger than 30mag. Evidence 

for hot dust from ISOCAM only. 

3C020: Evidence for hot dust by ISOCAM. 

3C022: Fitting the visible and ISOCAM data points by a 

Planck function requires a rest frame temperature of 3000 K. 

If one allows for moderate reddening of Ay = 1 mag in the 

galaxy the temperature would rise to 3600 K. Such a value is 

compatible to a normal stellar population. However, near IR 

imaging (Simpson et al. 1999) suggests a possible detection of 

dust. 

3C031: ISOCAM emission comes from a compact region of 
about 8" whereas the shown optical data are the magnitudes 
of the whole galaxy (180"). The MIR emission of the stars 
are therfore only a small contribution to the ISOCAM fluxes. 
Evidence for hot dust by ISOCAM. 

3C033.1, 3C048 and 3C061.1: Evidence for dust by ISOCAM 
and ISOPHOT. 

3C66B: Evidence for hot dust comes from ISOCAM data. 
3C079: Evidence for dust from ISOCAM and ISOPHOT. 
3C083.1: Star superposed just 3" southeast of the nucleus 
makes modelling of this galaxy impossible. 
3C098: Evidence for hot dust from ISOCAM only. 



3C231: M82 starburst model as described by Kriigel & 
Siebenmorgen (1994). 

3C249.1: Evidence for dust from ISOCAM and ISOPHOT. 

3C265: Evidence for dust from ISOCAM only. 

3C270: No evidence of dust from ISOCAM. 

3C272.1: Evidence for dust from ISOPHOT only (Haas et 

al. 2001). 

3C273: Possible hot dust detection by ISOCAM. 

3C274: No evidence of dust from ISOCAM. 

3C277.3: Evidence for dust from ISOCAM only. 

3C286: Evidence for hot dust by ISOCAM only. 

3C288.1: IRAS 60 /im maybe due to confusion. Likely detection 

of hot dust by ISOCAM. 

3C293: Evidence for hot dust from ISOCAM only. 
3C295: Evidence for dust from ISOCAM and ISOPHOT. 
3C296: IRAS 100 /im maybe due to confusion. Likely detection 
of hot dust by ISOCAM. 

3C303.1: Fitting the visible and ISOCAM data points by a 
Planck function requires a rest frame temperature of 2600 K. 
If one allows for moderate reddening of Ay = 1 mag the tem- 
perature would rise to 2900 K. Such a value is still lower than 
the expected mean stellar temperature in the galaxy. Evidence 
for dust from ISOCAM and ISOPHOT. 
3C305: Evidence for dust by ISOCAM and ISOPHOT. 
3C309.1: Evidence for dust by ISOCAM and ISOPHOT. 
3C319: Likely detection of hot dust by ISOCAM. 
3C321: Evidence for dust from ISOCAM and ISOPHOT. 
3C324: Evidence of hot dust from ISOCAM only. 
3C330: Evidence for dust from ISOCAM only. 
3C332: Evidence for hot dust from ISOCAM only. 
3C336: Evidence for hot dust from ISOCAM only. 
3C338: Evidence for dust from IRAS only. 
3C341: Evidence for dust from ISOCAM only. 
3C343: Possible detection of dust by ISOCAM only. 
3C345: Evidence for hot dust from ISOCAM only. 
3C34&. Evidence for dust from ISOCAM only. Large deviation 
with ISOPHOT 12 ^m fluxes given by Haas et al. (2001). 
3C349: Evidence for dust from ISOCAM only. 
3C351: Evidence for dust from ISOCAM and ISOPHOT. 
3C356: Evidence for dust from ISOCAM and ISOPHOT. 
3C371: No evidence for dust since the contribution of the syn- 
chrotron radiation to the ISOCAM data is unclear. 
3C379.1: Evidence for dust by ISOCAM only. 
3C380: The measured radio core flux is similar to the total ra- 
dio emission, so that the synchrotron component, interpolated 
to the IR, is strong. Still the ISOCAM detection gives evidence 
of dust emission. 

3C381: Evidence for dust from ISOCAM and ISOPHOT. 
3C386: Evidence for dust from ISOCAM at 14.3/im only. 
3C401: In the infrared/submillimeter there are only upper lim- 
its which hints to model this galaxy. 

3C402: ISOCAM emission comes from a compact region of 
about 9" whereas the shown optical data are the magnitudes of 
the whole galaxy (55"). The MIR emission of the stars is there- 
fore only a small contribution to the ISOCAM fluxes. Evidence 
for dust from ISOPHOT (Haas et al. 2001) only. 
3C410: Evidence for dust from ISOCAM only. 
3C442'- ISOCAM emission comes from a compact region (< 
5") whereas the shown optical data are the magnitudes of the 
whole galaxy (130"). The MIR emission of the stars is therefore 
only a small contribution to the ISOCAM fluxes. Weak 12/im 
detection might be still photospheric. No evidence for dust. 
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3C44& ISOCAM emission comes from a compact region of 
about 7" whereas the shown optical data are the magnitudes of 
the whole galaxy (90"). The MIR emission of the stars is there- 
fore only a small contribution to the ISOCAM fluxes. Evidence 
for dust comes from ISOCAM only. 
3C449: Evidence for hot dust from ISOCAM only. 
3C452: Evidence for dust comes from ISOCAM only. 
3C456: Fitting the visible and ISOCAM data points by a 
Planck function requires a rest frame temperature of 2300 K. 
If one allows for moderate reddening of Av = 1 mag the tem- 
perature would rise to 2600 K. Such a value is still lower than 
the expected mean stellar temperature in the galaxy. Evidence 
for dust from ISOCAM only. 

3C459: Evidence for dust from ISOCAM and ISOPHOT. 
3C465: ISOCAM emission comes from a compact region of 
about 10" whereas the shown optical data are the magnitudes 
of the whole galaxy (180"). The MIR emission of the stars 
is therefore only a small contribution to the ISOCAM fluxes. 
Evidence for hot dust from ISOCAM only. 
3C469.1: Evidence for dust from ISOCAM and ISOPHOT. 
4C+72.26: Our highest redshift quasar (z=3.532). The photo- 
spheric contribution is unknown, so that there is yet no evi- 
dence for dust. The model is not well constrained by the sparse 
infrared/submillimeter data. 
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Fig. B.l. ISOCAM contours overlayed on DSS grey scale images centered on NED position of the object marked in 
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Fig.B.l. - continued -: 3C079 at 12.0^m (LW10): ±89.6, ±106.7, 213.5, 320.2, 384.3. 3C083.1 at 4.5 /zm (LW1): 
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Fig.B.l. - continued -: 3C303.1 at 12.0/xm (LW10): ±5.0, ±10.0, 17.9, 117.5, 176.3 3C305.0 at 12.0^ (LW10): 
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Fig.B.l. - continued -: 3C465 at 6.7^m (LW2): ±9.2, ±9.9, 19.7, 29.6, 35.5. 3C469.1 at U.Ofim (LW10): ±10.3, 
±20.7, 37.2, 181.7, 272.5. 4C+72.26 at 12.0 (LW10): ±0.4, ±0.7, 1.3, 114.0, 171.1. NGC382 at 12.0/im 
(LW10): ±27.3, ±54.6, 98.3, 173.4, 260.1. PMNJ0214-11582 at 12.0^m (LW10): ±4.6, ±9.1, 16.4, 250.0, 375.0. 
PCd91212130.83+245139.2 at 6.7/im (LW2): ±48.2, ±72.3, 235.9, 353.9, 424.6. 
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Table B.l. ISOCAM observations of 3CR sources ordered by TDT number. 
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